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We have found an exact analytical solution of the Bogoliubov-de Gennes equations for the 
Tkachenko modes of the vortex lattice in the lowest Landau level (LLL) in the thermodynamic 
limit at any momenta and calculated their damping rates. At finite temperatures both Beliaev and 
Landau damping leads to momentum independent damping rates in the low-energy limit, which 
shows that at sufficiently low energies Tkachenko modes become strongly damped. We then found 
that the mean square fluctuations of the density grow logarithmically at large distances, which in- 
dicates that the state is ordered in the vortex lattice only on a finite (although exponentially large) 
distance scale and introduces a low-momentum cut-off. Using this circumstance we showed that at 
finite temperatures the one-body density matrix undergoes an exponential decay at large distances. 

I. INTRODUCTION 

The physics of rapidly rotating Bose-condensed gases is governed by a collective behavior of nucleated vortices 
(see [lH3| for review). When the rotation frequency f2 along the axis perpendicular to the plane of rotation (x,y) is 
close to the trapping frequency uj in the x and y directions, the harmonically trapped Bose gas becomes essentially 
two-dimensional (2D) and can be described as a system of interacting bosons in the lowest Landau level (LLL). If the 
number of vortices is much smaller than the number of particles, the system is in the so-called "mean- field" quantum 
Hall regime Q, where vortices arrange themselves in a lattice with an intervortex spacing of the order of the "magnetic 
length" / = h/y/mfl. Under an increase in f2 the vortex lattice should melt and strongly correlated Quantum Hall 
states should emerge. 

The vortex lattice in the LLL or close to this regime has been obtained experimentally and studied theoretically 
assuming the presence of a macroscopic wave function \l/(r) (see [H for review and refs. in @). Presently, a new 
generation of experiments with rapidly rotating quantum gases is being set up. It is based on the use of artificial 
gauge potentials, which are obtained by means of combinations of laser fields and mimic the rotation of the system 
(see [l| for review). They are supposed to provide the "rotation" with frequency that is extremely close to uj, thus 
allowing for the observation of the lattice melting and for the emergence of strongly correlated states. This puts 
forward the question of correlation properties in the vortex lattice regime and revives the interest to excitation modes 
in this regime. 

The theory of elastic oscillations of a vortex lattice in incompressible superfluids has been constructed by Tkachenko 
@, and a linear spectrum of the wave dispersion has been obtained. The effect of a finite compressibility has been 
discussed in Refs. CHI, and it has been shown that the compressibility changes dramatically the dispersion 
relation in the low-momentum limit. The dispersion becomes quadratic due to hybridization with sound waves. The 
first experimental observation of Tkachenko modes in harmonically trapped Bose-condensed gases has been reported 
in Ref.Q. Theoretical analysis of these modes has been done in the mean- field approach for both the geometry of an 
infinite plane (Q — uj) and in the presence of remaining trapping 

mis- 

The calculations provided an explanation 
of the mode frequencies observed in the experiment Q ■ Theoretical studies based on the microscopic effective action 
showed the absence of long-range order at T = in the thermodynamic limit (infinite plane geometry 1 ) 12 ill , and the 
hydrodynamic approach revealed an algebraic decay of the one-body density matrix at large distances [15j . 

In this paper we find an exact analytical solution of the Bogoliubov-de Gennes equations for the Tkachenko modes 
of the LLL vortex lattice in the thermodynamic limit at any momenta and calculate their damping rates. Importantly, 
at finite temperatures both Beliaev and Landau damping mechanisms lead to momentum independent damping rates 
in the low-energy limit. This means that for sufficiently low energies the Tkachenko modes become strongly damped, 
which is consistent with the experimental results [5(. Using the obtained results for the excitation wavefunctions we 
then calculate the mean square fluctuations of the density and the one-body denity matrix. The density fluctuations 
grow logarithmically at large distances, which indicates that the state is ordered in the vortex lattice only on a finite 
(although exponentially large) distance scale and introduces a low-momentum cut-off. Using this circumstance we 
show that at finite temperatures the one-body density matrix undergoes an exponential decay at large distances. 
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II. THE GROUND STATE AND EQUATIONS FOR THE EXCITATIONS IN THE LLL 

APPROXIMATION. 



We consider a zero-temperature two-dimensional (2D) system of bosonic atoms in a harmonic trapping potential 
V(r) — mco 2 r 2 /2, rotating with frequency Q, around the axis perpendicular to the (x,y) plane. In the rotating frame 
the Hamiltonian has the form: 



H = d 2 r 



a 2 

ft + IftftM + VMft-th - Q.ftt z ^ 
2m 2 



(1) 



where tp(r) is the field operator, p is the momentum operator, m is the atom mass, L z is the operator of the 
orbital angular momentum, and g is the coupling constant for short-range atom-atom interaction. The non-linear 
Schroedinger equation for i/)(r, t) reads: 

iK^ = ^-i> + g^i)i> + V{T)i>-SlL z i>. (2) 
at 2m 

It is commonly assumed that in the mean-field Quantum Hall regime all particles are in the same macroscopic 
quantum state described by the wavefunction vF JJQ. In the ground state this wavefunction has the form ^(r,^ = 
^o(r) exp(— ifit), where /i is the chemical potential. The wavefunction 5'o( r ) is then governed by the Gross-Pitaevskii 
equation: 

D 2 

o + </|* | 2 *o + V(r)* - fi£,tf = M^o, (3) 

2m 

where it is normalized to the total number of particles N. For Q close to w, the ground state of the system corresponds 
to the LLL and, hence, the macroscopic wave function x I'o( r ) takes the form: 

*o(r) = V^fo(z)e-^ /2 , (4) 

where n — N/S is the mean density, with S being the surface area, z,z — (x ± iy)/l, and the function fo(z) is 
analytical in the {x,y) plane, with fo(z) exp(— |z| 2 /2) normalized to unity. Equation (J3J is then solved numerically 
by expressing as a superposition of the LLL single-particle eigenfunctions z™exp(— \z\ 2 / 2) / im\ . Alternatively, 
one solves a projected Gross-Pitaevskii equation obtained by acting on Eq. <j3j) with the LLL projection operator P. 
This operator acts on an arbitrary function F(z, z) as 

PF(z,z) — — [ dwdwexpl— \w\ 2 + zw]F(w, w), (5) 

7T J 

and transformes it to the function in the LLL. 

For O = uj we have the geometry of an infinite plane. The LLL is infinitely degenerate and the projected Gross- 
Pitaevskii equation reads 



Ng 

IT 



dwdwe~ 2w *>+^\f (w)\ 2 f (w) = /*/„(«), (6) 



where fx = /j,— Ml. 

The ground state solution of Eq. ([6]) is a triangular vortex lattice and the function fo(z) is expressed through the 
Jacobi Theta- function $1 Q: 

f (z)^(2v) 1 / 4 MV^z,q)e z2 / 2 , (7) 

with q = exp(i7rr), r = u + iv, v — Vo/2, u = —1/2. The chemical potential is then equal to jx = ang, with 
a = \/«E m , P (-l) mp exp{-7n;(TO 2 + p 2 )} = 1.1596. 

An alternative mean-field approach does not apriori assume the presence of long-range order and is based on the 
density-phase formalism. It is commonly employed for 2D Bose-condensed gases at finite temperatures and for weakly 
interacting ID bosons, where the long-range order is destroyed by long-wave fluctuations of the phase. The key 
condition of this approach is related to small fluctuations of the density. One writes the field operator in the form 

Tp = exp i$ Vfi; = v^exp (— z$) (8) 
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with h and $ being the density and phase operators satisfying the the commutation relation [n(r), <p{v')\ = i5(r — r'). 
Representing these operators as n = no(r) + Sn and <f> = $o( r ) + one writes Eq. ([2]) in terms of the density and 
phase keeping only zero and first order terms in small fluctuations Sh and V<5$. To zero order we then have Eq. @ 
for *o( r ) = \/ n o{ r ) exp(i$o( r )) and, hence, the projected equation (5) for the function fo(z) introduced by Eq. (4). 
The solutions of two equations that are linear in 5n and V(5$ are obtained by representing these quantities in terms 
of elementary excitations characterized by the wavenumber k = {k x , k y }: 

Sh = -y/rioe - ' 2 ' ^ 2 ^[uk exp[— i<5> ] — u k exp[i$ ]] exp[— ie^t]^ + h.c. (9) 

k 



_j e -W 2 /2 „ 

(5$ = — - ) [uk exp[-i$ ] + w k exp[i$ ]] exp[-ie k £]ak + h.c, (10) 

k 

where a k , a k are operators of anihilation/creation of the excitations, and the functions iik> £k satisfy the Bogoliubov- 
de Gennes equations. In the LLL approximation they are analytical functions of z and follow from the projected 
Bogoliubov-de Gennes equations: 

2 3 P(|* |V)-.gP(*o«k) = (A + £kK 

2gP(\y \ 2 ik) - gP(*l<) = (A-ek)«k- (11) 



III. EXACT SOLUTION OF THE PROJECTED BOGOLIUBOV-DE GENNES EQUATIONS 

In this section we will use coordinates expressed in units of I, wavevectors in units of and energies in units of 
ng. The exact solution of Eqs. (fTTj) for the functions u^, reads: 



Clk 



u k = \z+-^±) e^-^ 2 e- k ' / 4 = ci k P(/oe lkr ) 



C 2 k , , 

Vk = 7$ fo{z ' 



ik 
~2 



+ \ e -ik-,/2 e -k*/* = C2k p( /oe - 1 



kr\ 



(12) 
(13) 



where k± = k x ±ik y , and the coefficients cik, C2k follow from the normalization condition f dx dy(\uk\ 2 — |i> k )e = 1. 
They are given by 



Clk 



£(k) + e k 
2e k 



1/2 



(14) 



C2k 



g(k) - ^ 
2e k 



1/2 



l^2(k)| fc2/8 

X- 2 (k) 



(15) 



with 



KQl) = 2KiQl)-K 1 (P), 



(16) 



nm — 7ri;(n 2 +m 2 ) 



C 



/7TU a; 77, + 2 ^7TU fe y 771 — k x / 4 



(17) 



*T 2 (k) = ^ £ (-l) mn e" 



(18) 



and 



7^(0) = JT 2 (0) = if(0) = a ~ 1.1596. 



(19) 
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The sums in Eqs. (|T71) and (IT8|) can be expressed in terms of Jacobi Theta-functions: 



Ki0c) ^ (^^) * tM^h^i (20) 



^ 2 (k) 



2u 



_g 2 (k y -\-ik x ) 



+1?2 



(fca; - iky ) ^77 

2> ' 



e 4 » 



(k x - iky)^ 



The spectrum of excitations has the form (see Fig. 1 and Fig. 2): 

e\ = |2K 1 (k)-^o| 2 -|^ 2 (k)| 2 , (21) 
It is consistent with the spectrum obtained in the hydrodynamic approach 0, [IH [111 , but we clearly identify the 





FIG. 1: Excitation energy e(k x , k y ) in units of ng. 



e 




k 



0.5 1.0 1.5 2.0 



FIG. 2: Anisotropy of the dispersion relation (|21[) . The upper curve is e(k,Q), the lower curve is e(k cos ^, fc sin ~), and the 
directions a; and y are chosen in the same way as in Fig. 1. 



anisotropy: for a given k the excitation energy is maximal in the x direction, and minimal at an angle of 30 degrees 



5 



from the x axis, (see Fig. 2). In the limit of small k, using the expansion of the functions Ki(k) and ^(k): 



Kx = a 



G4 



k 2 k A 

K 2 = a[l 1 

4 32 



where 



g E m , P (-l) mp mVexp[-^(m 2 +p 2 
V E m , P (-!) mp exp[-7rv(TO 2 + p 2 )} 

we find a symmetric spectrum. Restoring the dimensions it reads: 



(22) 



ng(kl) 2 ~ 0.2628 ng(kl) 2 



0.8219, 



(23) 



which exactly coincides with the hydrodynamic result of Ref. 19]. The obtained excitations are commonly termed 
as Tkachenko modes, although they have a quite different dispersion relation compared to elastic oscillations of the 
vortex lattice in incompressible superfluids obtained by Tkachenko 

For atoms (tightly) confined to the quasi2D geometry with frequency ujq, the coupling constant for the interatomic 
interaction is g — 2^/2tt t/ 2 a/ ml®, where a is the 3D scattering length, and lo = \J K/mujQ. In the case of 87 Rb rotating 
in the (x, y) plane with frequency Q ~ 100 Hz and confined in the perpendicular direction with ujq ~ 300 Hz, one 
has rig /Ml ~ 0.1 at the 2D density n ~ 3 x 10 8 cm~ 2 , which justifies the LLL approximation. Then, low-energy 
excitations have frequencies below 1 Hz. Note that for this example the quantity v = rrnl 2 representing the ratio of 
the number of particles to the number of vortices, i.e. the filling factor, is large and we are well in the mean-field 
regime. 



IV. DAMPING RATES OF THE EXCITATIONS 



Let us now calculate the damping of these excitations, which is caused by the interaction term of the Hamiltonian 
(1), containing a product of four field operators. For finding the damping rate it is sufficient to use a linearized form 
of the field operator, i.e. put ip = ^/ (1 + Sn/2n + iS<&). Using equations (4), © and (JTUJ) we then have: 



= ^/nfo(z) + ^[u k a k exp(-ie k t) - Wk« k exp(ie k £)^ exp(-i/it - \z\ 2 /2). (24) 

We first consider the Beliaev damping mechanism 22] in which a given excitation with wavevector p decays into 
two excitations with lower energies and wavevectors (k and q). The part of the interaction Hamiltonian that causes 
the Beliaev damping contains three operators of the excitations and has the form: 

V = gVn^2 J dxd y [/o w p u k w q + 2/oWpUqWk - /o^p^kUq - 2/ *M p w k z; q ] e~ 2|z| a q a k a p + h.c. (25) 

k,q 

At a finite temperature T we have to take into account thermal occupation of the states with momenta k and q and 
the presence of the reversed process in which excitations with momenta k and q recombine into the excitation with 
momentum p. Using the Fermi golden rule the damping rate is given by 

r p = Y E K k > ql^lP)!^ 1 + N * + N i) S (tr - £ k - e q ), (26) 

k,q 

with -/V ki q = [exp(e kj q/T) — 1] _1 being equilibrium occupation numbers for the excitations. The excitation energy 
thus acquires the imaginary part and becomes e p — iM p /2. 

For performing the calculations it is convenient to use the functions u k , u k in the form of Bloch waves: 

u k (r) exp(-|z| 2 /2) = Ci ^o(z + ik+l/2) exp( _ fc2 p /g) exp ( ikr / 2 ) ) ( 2 7) 

v N 



«k(r) exp(-|z| 2 /2) = c 2k ^° (z ^_+ l l 2 ) e ^^ k 2 l 2 /8) exp(-ikr/2). (28) 

v N 



6 



In the low-energy limit where pi <C 1, the transition matrix element is equal to 



a 5/2 I Nng (- ek 6q 



S y 8e p efce 9 [ if (?) if (p) J " H 

where <5 P ,k+q is the Kronecker symbol. After a straightforward algebra Eq. (|29|) is reduced to 



/i n>i \ a ng 2 (k 4 + q 4 - p 4 )l 

(k ' q|y|p> = ivf^v — — ^ — ^ (30) 



Equation ([25)1 then yields: 



and we obtain: 



^^-ngl 2 (p 2 - k 2 - (p 2 + k 2 - 2fcpcos( 



r *> = IT- h 2 f dxx 2 ^/T^ coth x 2 

07TT) H Jo \2l 



(31) 



At T = we immediately find: 

a 9 2.nn„9 2 



^=128^ ^0.011^. (32) 



For the ratio of the damping rate to the excitation energy we then obtain: 

hT p0 _ 0.13 



(33) 



Cp IS 

In the mean- field regime we should have v — irnl 2 3> 1, since this quantity (filling factor) represents the ratio of 
the number of particles to the number of vortices. Therefore, we have hT p o <C e p at any p. Thus, at T = Tkachenko 
modes are good elementary excitations in the entire mean-field Quantum Hall regime. 

The situation changes drastically at finite temperatures. Equation (f3T|) yields T p o + T p t, and for e p -C T the 
temperature-dependent part of the damping rate is independent of p and proves to be 

r ^=4^?^ Lo 4^ ^ <<T - (34) 

For excitation energies e p 3> T/v we have KT p t <C e p , i.e. the damping rate is small and, hence, Tkachenko modes 
are good elementary excitations. Moreover, for e p > T the damping rate starts to decrease with increasing p. In 
particular, for e p ^> T equation (|3ip gives: 

r pT = ^C(3/2)^-f-); e p »T, (35) 



4^3/2 fa 



where £(3/2) is the Riemann zeta- function, and the imaginary part of the excitation energy is negligible compared to 
the real part e p . However, excitations with energies 

T 

£p^ e c = — (36) 
v 

are overdamped. 

Note that our analysis was assuming the so-called collisionless regime, where pumping the mode with a given 
energy e p one does not disturb the equilibrium distribution function Nk for thermal excitations involved in the 
damping process. This means that the relaxation of their distribution function occurs on a time scale tr h/e p . For 
the discussed Beliaev damping, the thermal excitations that are involved in the damping of the mode with energy e p 
also have energies ~ e p . So, we have tr ~ T^ 1 , an d excitations with energies e p ^> e c are well in the collisionless 
regime. However, for e p < e c we have the condition tr£ p /H < 1, and these excitations enter the hydrodynamic regime 
(see [23|). A dimensional estimate for their damping rate is T ~ T p (e p Tji/h) : with T p being the damping rate in the 
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collisionless regime. Thus, for the modes with energies e p < e c we have the damping rate approaching e p , i.e. they 
are significantly damped. 

It should be noted that at finite temperatures we also have the Landau damping in which a given excitation with 
momentum p interacts with a thermal excitation (momentum k) , both are annihilated and an excitation with a higher 
energy and momentum (q) is created. There is a reversed process as well. The interaction Hamiltonian that causes 
this damping is: 

Vl = gVn^2 / dxdy [/o u q u k" P + 2/o{£6 q iip - fovlijpVq - 2/ it q 6*M k ] + h.c, 
and for the damping rate we have 

k,q 

The calculations are similar to those made above for the Beliaev damping. In the low-energy limit in both limiting 
cases, f p « T and e p ^> T, the results are the same as in the case of Beliaev damping, but with a twice as small 
numerical coefficient. Thermal excitations that are involved in the damping of the mode with energy e p also have 
energies ~ e p . We thus see that the Landau damping does not change our conclusion made from the analysis of the 
Beliaev damping. Namely, excitations with energies e p 3> e c are well in the collisionless regime, with the damping rate 
T p <C e/H. On the other hand, excitations with energies e p < e c enter the hydrodynamic regime and are significantly 
damped. 

For temperatures of the order of tens of nanokelvins and filling factors v of the order of hundreds, like in experiments 
, 0] where the LLL regime has been reached, Tkachenko modes with frequencies of the order of 1 Hz or lower should 
already be in the regime of strong damping. The strong damping of Tkachenko modes in this frequency range has 
been observed in the JILA experiment [5j, although this experiment was not yet in the LLL regime. 



V. ONE-BODY DENSITY MATRIX 



We now discuss correlation properties of rapidly rotating bosons in the mean-field Quantum Hall regime at zero 
and finite temperatures. For this purpose we will use the field operators in the form (JSj> . Due to small fluctuations of 
the density the one-body density matrix takes the form: 

9l (r) = (^t(r)^(O)) = *5to*o(0) exp {-~<(<5$(r) - <5$(0)) 2 ) j . (37) 

In the low-momentum limit where M « 1, omitting the term ik + l/2 in the argument of Wo in equations (1271) and 
(|28[). the operator of the phase fluctuations given by Eq. (flO|) becomes: 

6Q(r) = - % - (cik t_ C2k) exp(zkr/2) a k + h.c. (38) 
Then, using equations (|14[) and (1151) for the mean square fluctuations we obtain: 

((<5$(r) - ^(0)) 2 ) =*9 (1 + 2jVfc) [l ~ Mkr/2)], (39) 

J (2tt)^ e fe 

where Jo is the Bessel function. 

At T = using the low-energy spectrum l|23|) equation (|39l) immediately gives: 

o i rPo fjif 

<(**(r) - <5$(0)) 2 ) = -=- !T [1 - Mkr/2)}. (40) 

v J k 

The upper bound of the integration in Eq. (|41[) is pq ~ l/l, which represents the boundary of the first Brilluen zone, 
and for r 3> I we find 

(0$(r) - ^(0)) 2 ) * -j= \ In (^j , (41) 
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with C — 0.5772 being the Euler constant. For the density matrix we then have an algebraic decay at large distances: 

SiW« - , r»l, (42) 



which reproduces the result of Ref. 15]. So, in the thermodynamic limit there is no long-range order even at T = 0, 
and we are dealing with a phase-fluctuating Bose-condensed state. 

At finite temperatures we have to take into account thermal fluctuations of the phase. Before finding g\ (r) we 
calculate the mean square fluctuations of the density ((Sn(r) — 5h(0)) 2 )/n 2 , which should be small for the validly of 
the mean-field approach. This is the case at T = 0, but at finite temperatures thermal fluctuations drastically change 
the situation. Using equations ((9J, (|27| . and ([28]) in the low-momentum limit we have: 

<(fft(r)-W(0))»> _ f d>k ga#P Jo(W2)](1+2jVfe) . (43) 



o<fc<i-i ( 27 2 2e fc 

For obtaining this relation we put c\k — C2k — ct and expanded ^q(z + ik + 1/2) in powers of k up to the first order. 
Then, omitting small vacuum fluctuations and writing Nk ~ T/e^ we reduce Eq. (|43[) to 

« a '' (r) -f"" )2 > - f ' [1 - Jo (W2)]f , (44) 

(xqngv J k 

where / = Z for T 3> w?, and for T ng we have / = k^ 1 with the momentum fcy following from the condition 
6k T = T. The integration is straightforward and it yields: 

((5n(r) - 6n(0)) 2 ) 4T ^ /_r\ 



n a.T]ngv \2l 

i.e. the density fluctuations grow logarithmically with the distance. This means that at finite temperatures the 
mean-field approach can be employed only on a distance scale r < r$ = 2lex^{arjngv/'iT). Then the mean square 
fluctuations of the density are small. In other words, the state is ordered in the vortex lattice only at r < ro- 

This introduces a low- momentum cut-off r^" 1 . For the thermal mean square fluctuations of the phase we then have: 

7T Jra 1 e k exp(e fc /i ) - 1 

The most important contribution to the integral comes from low momenta, so that we can expand the exponent in 
the denominator of Eq. (|46[) and put the upper limit of integration equal to infinity. Then Eq. (|46|) takes the form: 



<(<H>(r) - «$<& (0)) 2 ) = / HI [1 _ J (fcr/2)], (47) 

mjngv J r -i k 6 

A straightforward integration at distances r 3> I yields: 

((<5$(r) - 5mf) = -^-^ b . (48) 
Omitting vacuum phase fluctuations we then obtain an exponential decay of the density matrix: 



gi(r) tx exp 



4T r\ (r (T) 
■ m 



arjngv I 2 



(49) 



For systems with a finite size L < ro we have to replace ro with L in Eq. (|49[) . 

Stricktly speaking, equation (|4"9"|) is applicable only at distances r <C r , where the system is ordered in the lattice. 
However, we clearly see that for r approaching ro the density matrix practically drops to zero and, hence, it should 
remain close to zero at larger distances. 
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VI. CONCLUDING REMARKS 

Concluding our work we would like to make a few remarks. First of all, the length scale ro on which the system 
is ordered in the vortex lattice at finite temperatures is exponentially large and in cold atom experiments it exceeds 
the size of the sample. Even on approach to the melting point, assuming T ~ ng and v ~ 20 we have r ~ 3002. 
However, the effect of finite temperature is likely to be dramatic for the visibility of the vortices. Equation (1451) shows 
that even fairly well in the mean-field regime, for example at v ~ 40 the mean square fluctuations of the density are 
((Sh(r) — 5n(0)) 2 ) ~ 0.3 for r ~ 101 and T ~ ng. This should significantly reduce the visibility of the vortices. We do 
not claim that these arguments explain the reduction of the vortex visibility in the ENS experiment Q, but rather 
attract attention to this finite-temperature effect for future studies. 

It is also worth mentioning that at finite temperatures the damping of Tkachenko modes may serve as a signature 
of the approach to the melting point of the lattice. The characteristic excitation energy e c ~ T/v below which these 
modes are strongly damped, increases significantly with decreasing the filling factor and becomes of the order of 10 
Hz for v — 40 even at temperatures as low as 20 nK. 
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